Published online 26 March 2013 



Nucleic Acids Research, 2013, Vol. 41, No. 9 5075-5089 

doi:10.1093/nar/gkt203 



Structural complexity of Dengue virus untranslated 
regions: c/s-acting RNA motifs and pseudoknot 
interactions modulating functionality of the 
viral genome 

Joanna Sztuba-Solinska 1 , Tadahisa Teramoto 2 , Jason W. Rausch 1 , Bruce A. Shapiro 3 , 
Radhakrishnan Padmanabhan 2 and Stuart F. J. Le Grice 1 * 

1 RT Biochemistry Section, HIV Drug Resistance Program, Frederick National Laboratory for Cancer Research, 
Frederick, MD 21702, USA, department of Microbiology and Immunology, Georgetown University School of 
Medicine, Washington, DC 20057, USA and 3 CCR Nanobiology Program, Computational RNA Structure Group, 
Frederick National Laboratory for Cancer Research, Frederick, MD 21702, USA 

Received January 16, 2013; Revised February 11, 2013; Accepted March 5, 2013 



ABSTRACT 

The Dengue virus (DENV) genome contains multiple 
c/s-acting elements required for translation and rep- 
lication. Previous studies indicated that a 719-nt 
subgenomic minigenome (DENV-MINI) is an efficient 
template for translation and (-) strand RNA synthe- 
sis in vitro. We performed a detailed structural 
analysis of DENV-MINI RNA, combining chemical 
acylation techniques, Pb 2+ ion-induced hydrolysis 
and site-directed mutagenesis. Our results highlight 
protein-independent 5-3 terminal interactions 
involving hybridization between recognized c/s- 
acting motifs. Probing analyses identified tandem 
dumbbell structures (DBs) within the 3' terminus 
spaced by single-stranded regions, internal loops 
and hairpins with embedded GNRA-like motifs. 
Analysis of conserved motifs and top loops (TLs) 
of these dumbbells, and their proposed interactions 
with downstream pseudoknot (PK) regions, pre- 
dicted an H-type pseudoknot involving TL1 of the 
5' DB and the complementary region, PK2. As 
disrupting the TL1/PK2 interaction, via 'flipping' 
mutations of PK2, previously attenuated DENV 
replication, this pseudoknot may participate in regu- 
lation of RNA synthesis. Computer modeling implied 
that this motif might function as autonomous 
structural/regulatory element. In addition, our 
studies targeting elements of the 3' DB and its 
complementary region PK1 indicated that 



communication between 5-3 terminal regions 
strongly depends on structure and sequence com- 
position of the 5' cyclization region. 



INTRODUCTION 

Unlike cellular mRNAs, which only contain translational 
regulatory motifs, viral RNA genomes have evolved 
functional regions, which act at different stages of their 
life cycle. These specific structures, their interactions 
and association with host proteins and ligands, together 
orchestrate the multifunctionality of viral RNA genomes. 
Often, RNA-RNA interactions spanning significant dis- 
tances (1,2) bring regulatory motifs into proximity, 
changing our view of a functional viral RNA genome 
from that of a linear molecule to one whose 3D structure 
is an important contributor during the life cycle (3,4). 

Organization of the Dengue virus (DENV) genome 
represents an example of intra-genomic long distance 
interactions that modulate molecular processes. DENV 
is a mosquito-borne, single-stranded, positive (+)-sense 
RNA virus belonging to the Flaviviridae family, whose 
members are responsible for diseases of major health 
importance, such as yellow fever, Japanese encephalitis, 
tick-borne encephalitis, dengue fever and a dengue hem- 
orrhagic fever. The DENV genome encodes a single-open 
reading frame, flanked by highly structured untranslated 
regions (UTRs) (5-7). The secondary structure of DENV 
UTRs has attracted attention, as they have been shown to 
encompass motifs involved in regulation of translation, 
replication, transcription and viral pathogenesis (7-10). 
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It has been proposed that the 3'-localized elements, 
including the terminal stem-loop (3' SL) together with 
cyclization sequence (3' CS), the upstream of AUG 
region (3' UAR) and the downstream of AUG region (3' 
DAR), are crucial for the initiation of (— ) strand, RNA- 
dependent RNA synthesis (2,5). These motifs are 
proposed to interact with their complementary elements 
positioned at the 5' terminus, the 5' CS, the 5' UAR and 
the 5' DAR, resulting in long-range RNA RNA 
mediated genome circularization (5,6,11). Following hy- 
bridization of the 5' and 3' termini, the DENV genome 
has been proposed to endure extensive rearrangements of 
both UTRs, opening of the 3' SL to form an extended 
duplex with the 5' UAR and formation of a double- 
stranded region via 5' and 3' CS interactions (12). As a 
result, the viral NS5 polymerase initially docked at the 5' 
stem-loop A (SLA) would be relocated to the structurally 
reorganized 3' terminus, allowing launch of (— ) strand 
synthesis (6,11). Additional structural motifs of the 
3'-UTR have also been shown to modulate DENV molecu- 
lar processes. Phylogenetic comparisons have suggested 
that the tandem, almost identical dumbbell structures, 
designated 3' and 5' DB, might participate in regulation 
of RNA synthesis and translation (13-15). The leftmost 
hairpins of both DBs contain a set of identical five nucleo- 
tide motifs, referred to as TL1 and TL2 (5'-GCUGU-3'). 
The top loop (TL) regions have been suggested to pair with 
complementary motifs (PK2 and PK1, respectively) nested 
within A-rich sequences flanking both DBs (13,16). The 
TL1/PK2 and TL2/PK1 interactions would promote for- 
mation of two pseudoknots, which have been proposed to 
modulate both viral replication (15) and translation (8). 

Current knowledge regarding the secondary structure of 
DENV 3'-UTR derives from folding algorithms and 
phylogenetic comparisons (8,13). With respect to the 
structure of 5' terminal motifs, the majority of data ori- 
ginate from chemical or enzymatic probing of short RNA 
segments (12), often in the absence of their complementary 
partners (17). In this study, we have analyzed the interplay 
of these regions by solving the secondary structure of 
DENV 5' and 3' terminal regions in the context of a 
719-nt subgenomic minigenome (hereafter referred as 
'DENV-MINF) (18). We combined selective 2'-hydroxyl 
acylation analyzed by primer extension (SHAPE) with 
Pb 2+ ion-induced hydrolysis and site-directed mutagenesis 
of DEN V-MINI mutant RNAs to examine the previously 
proposed cw-acting RNA motifs and their potential 
binding partners, as well as identify long-range tertiary 
interactions between them. Our findings provide new 
insights into the secondary and tertiary interactions of 
DENV-MINI RNA and contribute to the understanding 
of viral RNA replication. 

MATERIALS AND METHODS 

DENV-MINI mutants 

Construction of the plasmid encoding DENV-MINI RNA 
has been described previously (18). Mutants were gen- 
erated by site-directed mutagenesis using the Agilent 
QuikChange Lightning Site-Directed Mutagenesis kit 



according to the manufacturer's protocol 
(Supplementary Table SI). 

In vitro RNA synthesis 

DNA templates for in vitro transcription were generated 
by PCR amplification of the plasmid DENV-MINI and 
the corresponding mutants, using primers listed in 
Supplementary Table SI. All PCR experiments were per- 
formed using Invitrogen Platinum® Taq DNA polymerase 
High Fidelity. Transcripts were synthesized with the 
Ambion T7-MEGAscript following the manufacturer's 
protocol, and RNAs were purified from denaturating 8 
M urea/5% polyacrylamide gels, followed by elution and 
ethanol precipitation. Purified RNAs were dissolved in 
sterile water and stored at — 20°C. 

Selective 2'-hydroxyl acylation analyzed by primer 
extension (SHAPE) 

Twenty picomoles of RNA was heated at 90°C for 3 min 
and slowly cooled to 4°C. The volume was adjusted to 
150 nl in a final buffer of 50 mM Tris-HCl (pH 8.0), 
100 mM NaCl and 5mM MgCl 2 . Samples were incubated 
at 37°C for 15 min. Folded RNA was divided into two 
equal portions (72 ul) treated with 8 ul of 30 mM 
/V-methylisatoic anhydride (NMIA) in anhydrous 
dimethyl sulfoxide (DMSO) (+) or DMSO alone (-). 
Tubes were incubated at 37°C for 50 min, and RNA was 
precipitated at — 20°C with 60ng/ul glycogen, 0.3 M 
sodium acetate (pH 5.2) and 3 volumes of cold ethanol. 
Precipitated RNA was collected by centrifugation, 
washed once in 70% ethanol and resuspended in 10 ul of 
water. Five picomoles of Cy5-labeled (for NMIA-modified 
samples) or Cy5.5-labeled (for unmodified samples) primer 
S720 or S440 (Supplementary Table SI) in 7ul H 2 0 was 
annealed to the RNA at 85°C for 1 min, 60° C for 5 min and 
35°C for 5 min. RNA was reverse transcribed at 50°C for 
20 min with 100 U reverse transcriptase (RT) (Invitrogen 
superscript III), lx RT buffer (Invitrogen), 5mM 
dithiothreitol (DTT) and 500 mM dNTPs (Promega). 
RNA was hydrolyzed with 200 mM NaOH (4 M) for 
5 min at 95° C, and reactions were neutralized with an 
equivalent volume of HC1 (2 M). Sequencing ladders 
were prepared using the Epicentre cycle sequencing kit ac- 
cording to the manufacturer's instructions and primers 
labeled with WellRed D2 and LicorIR-800 dyes. 
Modified and control samples were mixed with the 
sequencing ladders, precipitated as described earlier in 
the text, dried and resuspended in 40 ml of deionized 
formamide. Primer extension products were analyzed on 
a Beckman CEQ8000 Genetic Analysis Systems as 
described previously (19). Electropherograms were pro- 
cessed using the SHAPEfinder program, following the 
software developer's protocol and included the required 
pre-calibration for matrixing and mobility shift for each 
set of primers as described (20,21). Briefly, the area under 
each negative peak was subtracted from that of the corres- 
ponding positive peak. The resulting peak area difference 
at each nucleotide position was then divided by the average 
of the highest 8% of peak area differences, calculated after 
discounting any results greater than the third quartile plus 
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1.5x the interquartile range. Normalized intensities were 
introduced into RNAstructure version 5.3 (22). 

Pb 2+ ion-induced cleavage 

Before Pb 2+ ion-induced cleavage, RNA transcripts 
(lOpmol) were subjected to refolding as described earlier 
in the text. A freshly prepared solution of Pb(OAc) 2 was 
added to a final concentration of 0.5mM, and the reaction 
mixtures were incubated for lOmin at room temperature 
(23). Reactions were quenched by adding an equal volume 
of 20 mM EDTA followed by ethanol precipitation. 
Cleavage sites were detected by primer extension as 
described earlier in the text (19). 

Antisense interfered SHAPE 

Locked nucleic acid (LNA)/DNA chimeras were ordered 
from Exiqon, the sequences of which are provided in 
Supplementary Table SI. Chimeras were added at 
10-fold excess after folding. Samples were subsequently 
incubated at 37°C for 15min before NMIA treatment 
(see earlier in the text). To quantify alterations induced 
by ai-oligonucleotides, raw data were processed as 
described earlier in the text. Tertiary interactions were 
introduced manually into RNAstructure version 5.3 
generated models (22). 

3D modeling of DENV-MINI RNA 

The secondary structure of DENV-MINI RNA predicted 
by RNAstructure software version 5.3 (22) and chemical 
probing data from SHAPE (21) were used to generate 3D 
models for 474-nt DENV-MINI RNA by RNAComposer, 
version 1.0 (http://rnacomposer.cs.put.poznan.pl/). The 
quality of predicted models was evaluated using 
MolProbity tool (24). Where required, spatial positioning 
of RNA substructures were adjusted and the affected 
regions subjected to energy minimization using Discovery 
Studio v3.5 (Accelrys, San Diego, CA, USA). 

RESULTS 

Probing analysis of DENV-MINI RNA verifies the 
5'-3'-UTRs long-range interactions and suggests a TL1/ 
PK2 pseudoknot 

SHAPE has been developed to study RNA secondary 
structure by examining backbone flexibility (directly 
related to the base pairing) at each nucleotide position 
via reactivity with a specific electrophilic reagent (21). 
We applied this technique to the in v/fro-synthesized 
719-nt subgenomic 'minigenome' RNA (DENV-MINI) 
template, which contains the 5'-terminal 226 nt, 42 nt 
from the C-terminal coding region of NS5, including the 
UAG termination codon, and the 451-nt 3' terminus (18). 
Previous studies showed that this RNA molecule harbored 
essential ra-acting elements for efficient translation 
(25,26), as well as (— ) strand RNA synthesis in vitro by 
either DENV2-infected cell lysates or purified RNA- 
dependent RNA polymerase (RdRP) (27,28). DENV- 
MINI RNA was also used for secondary structure 
prediction analysis using the MPGAfold algorithm (29) 



to infer the structure-function role of the TLs and PKs 
regions in DENV replication and translation (8). 

Reactivity of NMIA at each nucleotide position of 
DENV-MINI RNA is shown on Figure 1A and B and 
Supplementary Figure SI. The most reactive and thus 
least conformationally constrained residues have a reactiv- 
ity >0.8 U and are depicted in red. Nucleotide positions 
with reactivity <0.2 U, indicative of fully base paired 
residues, are in black, whereas residues of intermediate 
reactivity are color-coded according to the key. Data 
shown are an average of at least three experiments. 

Minimal free-energy modeling using SHAPE data as 
pseudo free-energy constraints is consistent with 
proposed 5'-UTR motifs in the context of DENV-MINI 
RNA (Figure IB and Supplementary Figure SI) (6,17,30). 
Two regions that exhibit high sensitivity to NMIA modi- 
fications were located within the 5' stem-loop A (SLA), 
A31-U34 and U48-A50, which corresponded to the TL 
and the side stem-loop (SSL), respectively. Three helical 
regions of 5' SLA, U3-G9:C64-G70, CI 1-G15A.56-G61 
and C24-A29:U35-G40, were insensitive to NMIA. The 
remainder of the reactive nucleotides were situated within 
the U62-U63 bulge and internal loops G17-A19 and 
G22-A23. Downstream of SLA, the oligoU tract 
spanning nt U71-U80 was predicted to be single- 
stranded. This region was followed by the unreactive 
double-stranded upstream of AUG region (dsUAR), 
formed via a long-range interaction between the 5'- and 
3'-UTRs (nt A81-G96 and nt C638-U654). Here, the only 
reactive residues occupied the vicinity of bulge (C646), 
mismatch (G84) and the closing base pair of the helix. 
The internal loop enclosing the first AUG codon (nt 
A97-A104 and G631-C637) was predicted to be single- 
stranded followed by the double-stranded 6-nt long down- 
stream of UAR region (DAR) containing reactive residues 
within the mismatch (A 108) and the closing pairs of the 
helices. The capsid hairpin (cHP, A115-U133) (12), dis- 
played an unreactive stem followed by a 7-nt reactive 
apical loop. As suggested by previous chemical probing 
studies, the 5—3' terminal interactions were mediated by 
base pairing between the 5' and 3' cyclization sequences, 
which formed 11 -nt double-stranded cyclization region 
(CS) (U134-G144 and C614-A624) (12). Subsequently, 
two short hairpins that formed downstream of the 5'-3' 
CS were found within the sequence of the 5'-end of the 
capsid protein gene. One contains a purine-rich and mod- 
erately reactive apical loop (A152-U159) and the other 
with a reactive 3-nt loop and A-C mismatch (G165- 
C185). Both hairpins were followed by 5-nt helical 
region (A190-U194 and A444-U448), which provided 
additional linkage between the 5' and 3' termini. 

With regard to the 3'-UTR, SHAPE provided informa- 
tion on the 3' SL structure, obscured only at the extreme 3' 
terminus because of the presence of hybridized primer 
(Figure IB and Supplementary Figure SI). The 3' SL 
displayed NMIA sensitivity at the apical loop (A674- 
C680) and insensitivity at its leftmost side. The fact that 
this large 79-nt stem-loop was followed with reactive 
stretch of 3 nt (C655-U657) and unreactive dsUAR 
region (nucleotides A81-G96 and nucleotides C638- 
U654), enforced the proposed conformation. Previous 
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Figure 1. Secondary structure of the nucleotides 1-719 of DENV-MINI RNA. (A) Processed SHAPE reactivities of DENV-MINI RNA are pre- 
sented as a function of nucleotide position. Red and orange notations are expected to fall in single-stranded regions, whereas bases indicated in black 
and green correspond predominantly to either double-stranded regions or putative tertiary interactions. (B) Three distinct DENV RNA domains 
have been annotated, namely, 5'-3'-UTRs, DBs and VR. Nucleotide positions are numbered every 20 nt. The abbreviations correspond to the 
following as-acting RNA elements: UTR, untranslated region; SLA, stem-loop A; dsUAR, double-stranded upstream of AUG region; cHP, capsid 
hairpin; CS, cyclization sequences; DB, dumbbell; TL, top loop; SL, stem-loop. The pattern of Pb 2+ -induced hydrolysis obtained for interacting 
regions is indicated in the upper left insert. 
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phylogenetic comparative studies suggested that the 
3'-UTR enclosed two almost identical DBs that might 
assume a pseudoknot conformation with their down- 
stream complementary regions (13,16). SHAPE data con- 
firmed the presence of both, and in addition provided 
detailed insight into their secondary structure. The 3' 
DB contained a 6-bp unreactive stem (U536-G541 and 
C598-A603) and two hairpins on opposing ends (G547- 
U571 and G584-C597). The apical loop of the leftmost 
hairpin incorporated the highly reactive TL2 sequence 
(G558-U562). The reactivity of TL2, combined with the 
fact that its downstream complementary partner PK1 
(A613-C617) was involved in formation of 5'-3' CS 
region, would likely prevent the previously suggested 
TL2/PK1 pseudoknot interaction (13,16). Regarding the 
5' DB, it likewise contained an unreactive stem (G449- 
G454 and C511-C516) and two hairpins on opposing 
ends (C462-G482 and G497-C510). Here, however, the 
apical loop of the leftmost hairpin incorporating the 
TL1 sequence (G467-G477), although predicted to be 
single-stranded, harbored mainly unreactive nucleotides. 
The fact that its complementary partner PK2 was 
likewise unreactive and enclosed within single-stranded 
loop (G526-C530) was strongly suggestive of a TL1/ 
PK2 base pairing interaction (Figure 2) (13,16). In 
addition, both DBs comprised reactive bulges (A491- 
C496 and A578-A583) and moderately NMIA-sensitive 
purine-rich internal loops. One of the internal loops of 
the 3' DB (G572-A575) was insensitive to modification 
as it contained a -GGAA- motif. Such -GNRA- tetraloops 
represent a unique fold, which include a trans sugar edge/ 
Hoogsteen edge G"A pair between the first and the last 



residues, a hydrogen bond between the 2'-OH of G and 
the N7 of A and 3' stacking of all tetraloop bases except 
the G. This architecture would account for lack of reactiv- 
ity towards NMIA (31). 

The region upstream of the DBs is often referred to as 
the variable region (VR), as its sequence differs between 
DENV serotypes (Figure IB and Supplementary Figure 
SI) (33). Our probing analysis indicated its propensity to 
form long single-stranded stretches (A274-A291, U364- 
A376, G418-C424 and U439-A443). One of the loops 
(U439-A443), composed almost exclusively of adenine 
residues, was mainly NMIA unreactive. It was shown pre- 
viously that consecutive adenine bases in single-stranded 
regions are often involved in stacking interactions, which 
would account for their inaccessibility towards NMIA 
modification (34,35). The 3'-UTR region also contained 
two stable tetraloops (C352-G363 and C425-G438). As 
expected, the -GNRA- tetraloop (C425-G438), restricted 
by extensive stacking and hydrogen-bonding interactions, 
was NMIA-insensitive (31). On the other hand, the - 
CAAG- tetraloop (C352-G363) behaved differently. 
Here, only its GC residues showed little reactivity, which 
could be interpreted as formation of an intraloop G*C 
base pair, whereas its AA residues were reactive. The 
apical loops of two other hairpins (G377-C392 and 
G396-C412) showed strong NMIA reactivity. Also, an 
extended hairpin containing internal asymmetric loop 
(A296-U301 and U345-G346), a bulge (C307-A313) 
and A-rich, moderately reactive apical loop (A323- 
A331) was predicted to form between nucleotides G292- 
C350. Another long hairpin (U206-A267) with a reactive 
apical loop (A233-A241) and two reactive internal loops 




Figure 2. TL1/PK2 pseudoknot interaction in DENV-MINI RNA. (A) SHAPE-predicted hairpin-type pseudoknot interaction comprises two helical 
regions, including base paring between PK2 (G526-C530) and TL1 residues (G470-U474), as well as the stem of the leftmost hairpin (C462-G466 
and C478-G482), and three single-stranded loops: A468-G470, U474-C476 and C516-C525. (B) The 2D map of TL1/PK2 interaction forming H- 
type pseudoknot generated by PseudoViewer3 (32). 
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(G218-A222 and C252-A256; A227-A229 and G245- 
A247) enclosed a linker sequence that stretched from 
A227-G268. 

Our secondary structure determination of DENV- 
MINI RNA was furthermore validated by single-strand- 
specific Pb 2+ ion-induced hydrolysis (23) (Supplementary 
Figure S2). The observed cleavage pattern correlated with 
SHAPE-predicted single-stranded regions. The 5'-UTR 
regions highly sensitive to Pb 2+ cleavage included the 
vicinity of SLA internal loop (C14-G15), SLA SSL 
(U48-A49) and the polyU tract (U71-U80). The SLA 
TL (C30-U34), as well as the apical loop of cHP (A121- 
A127), was hydrolyzed weakly, agreeing with previous 
studies showing that in most cases, bulges and internal 
loops were more readily cleaved with Pb 2+ ions than 
apical loops (23). The 5'-3' CS (U134-G144) and 
dsUAR (A88-G96) were resistant to cleavage, as they 
were predicted to be double-stranded. In addition, the 
two hairpins formed within the sequence of the 5'-end of 
the capsid protein gene (A152-G153 and A174-C175) 
were moderately sensitive to hydrolysis, whereas the 
following single-stranded region was hydrolyzed readily 
(C199-A205). 

Pb 2+ ions moderately cleaved some of single-stranded 
stretches of the 3'-UTR (A277-C289), whereas others, 
particularly the A-rich ones, had isolated cleavage points 
(A369 and C375) or were resistant to hydrolysis (U439- 
A443) (Supplementary Figure S2). The two tetraloops of 
the 3'-UTR (C352-G363 and C425-G438) were resistant 
to hydrolysis, consistent with previous studies, indicating 
that tetraloops showed little or no Pb 2+ -induced hydroly- 
sis because of extensive stacking interactions (36). The 
apices of other hairpins (G325-A327, A384-C385 and 
G402-U408) were moderately reactive. The most 
pronounced Pb 2+ -induced hydrolysis occurred within the 
apex of the 3' SL (C673-C675), the internal loops of DBs 
(A491-G495 and C579-G582) and the loops of their right- 
most hairpins (U500-G505 and A588-G592). Lack of 
Pb 2+ cleavage within TL1 and its putative complementary 
partner PK2 was again suggestive of their base pairing 
interaction (Figure IB insert and Supplementary Figure 
S2). In contrast, TL2 (G555-U562) was hydrolyzed 
readily, whereas its PK1 partner, predicted to be 
involved in the double-stranded 5—3' CS interaction, was 
resistant to cleavage. 

Antisense-interfered SHAPE data support a novel TL1/ 
PK2 pseudoknot interaction 

To examine the TL1/PK2 tertiary interaction (13,16), anti- 
sense-interfered SHAPE (aiSHAPE) strategy has been 
applied (37). In this method, displacement of one strand 
of an RNA duplex by hybridizing an antisense oligo- 
nucleotide would disrupt long-distance interactions and 
would be characterized by enhanced NMIA sensitivity 
of the displaced nucleotides (Figure 3A). 

Two chimeric LNA/DNA heptanucleotides were 
hybridized to either TL1 or PK2 regions to disrupt the 
potential base pairing interaction with their respective 
partner (Supplementary Table SI). LNA/DNA hepta- 
nucleotide 1A was complementary to the single-stranded 



region A523-G529, including PK2, whereas LNA/DNA 
heptanucleotide IB was complementary to the loop of the 
leftmost hairpin of the 5' DB, C471-G477 and included 
TL1 region. 

Hybridization of LNA/DNA 1A enhanced reactivity 
within TL1 residues (Figure 3B and Supplementary 
Figure S3). Median reactivities for nucleotides A468- 
G473 were 0.27, whereas the corresponding values for 
DENV-MINI were 0.05. In addition, LNA/DNA 1A 
hybridization formed an extensive barrier to reverse tran- 
scription (A503-G529), supporting the specificity of our 
aiSHAPE strategy (Supplementary Figure S3). 
Hybridization of LNA/DNA IB increased the median 
reactivities of the PK2 residues G526-C530 to 0.38, 
whereas the corresponding values for DENV-MINI were 
0.05 (Figure 3C and Supplementary Figure S3). Here, 
again, the specificity of LNA/DNA IB hybridization to 
TL1 was substantiated by formation of an extensive stop 
during reverse transcription (C463-G477) (Supplementary 
Figure S3). In both cases, NMIA reactivity profiles for the 
entire 719-nt DENV-MINI indicated minor 'off site' 
changes mainly occurring within the DB structures, 
whereas the 5' portion of the DENV-MINI was not sig- 
nificantly affected (Supplementary Figure S3). Despite 
this, the RNAstructure algorithm, when provided with 
pseudo energy constrains retrieved from LNA experi- 
ments, predicted the identical DENV-MINI secondary 
structure (data not shown). However, we cannot exclude 
the possibility that global DENV-MINI RNA folding 
might be altered by LNAs hybridization on a higher 
level that chemical probing cannot address. Despite 
these minor discrepancies, aiSHAPE data lend further 
credence to TL1/PK2 pseudoknot (Figure 2). 

Mutations within PK2 region disrupt TL1/PK2 
pseudoknot 

Previously, Manzano et al. (8) suggested putative roles of 
the TL/PK interactions in DENV replication. To disrupt 
potential TL1/PK2 and TL2/PK1 base pairing, they 
introduced point mutations by 'flipping' the correspond- 
ing sequences (TLs 5'-GCUGU-3'^5'-UGUCG-3' and 
PKs 5'-GCAGC-3'^5'-CGACG-3) in the context of the 
DENV2 replicon expressing Rluc. This strategy showed 
that the PK2 F i ip mutation attenuated DENV Rluc repli- 
cation, but not translation, suggesting that PK2 region is 
crucial for viral RNA synthesis. As the DENV2 Rluc 
replicon encompassed the same cw-acting regulatory se- 
quences as the DENV-MINI RNA used in our studies, 
we next sought to address the structure-function relation 
of PK2 mutations in the context of DENV-MINI RNA. 

SHAPE analysis of PK2 Flip DENV-MINI mutant RNA 
indicated major reactivity changes within the 5' DB struc- 
ture (Figure 4 and Supplementary Figures S4 and S5). The 
apical loop of the TL1 -containing hairpin was rendered 
NMIA-sensitive with median reactivities for G470-U474 
of 0.55, whereas the corresponding values for WT DENV- 
MINI were 0.05. The apex of the rightmost hairpin 
exhibited lower reactivity for residues G502-A506. In 
addition, the 'flipping' mutation of PK2 increased its sus- 
ceptibility to acylation (C526-G530) validating disrupted 
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Figure 3. Structural responses of the TL1/PK2 pseudoknot motif to antisense oligomers. (A) aiSHAPE principal. Hybridization of an interfering 
oligonucleotide (green) to one partner of the proposed RNA duplex increases acylation sensitivity of its base-paired counterpart. (B) Influence of 
LNA/DNA oligonucleotide (orange) hybridized to PK2 on chemical reactivity of TL1 residues. (C) The effect of LNA IB hybridization to TL1 loop 
on acylation of PK2 nucleotides. The gray residues represent the formation of an extensive stop during reverse transcription at the place of LNA 
hybridization. The native traces (red) are compared with antisense-interfered DENV-MINI RNA (violet). Nucleotide positions exhibiting increased 
NMIA reactivity in the presence of given LNA are distinguished by orange brackets. 



TL1/PK2 base pairing. Other regions did not exhibit sec- 
ondary structure rearrangements. Thus, data of Figure 4 
and all subsequent illustrations representing RNA struc- 
tures of DENV-MINI mutants present only the affected 
RNA motifs. 

Single-strand-specific Pb 2+ -induced cleavage of the 
PK2 F | ip construct verified that the 'flipping' mutation 
released both TL1 and PK2 from a base pairing inter- 
action, as both regions were now strongly susceptible to 
hydrolysis (Figure 4 insert and Supplementary Figure S6). 



In addition, the apex of the rightmost hairpin showed 
lowered cleavage intensity as compared with DENV- 
MINI RNA (Supplementary Figure S2). 

Mutations within PK1 region impair 5'-3' terminal long- 
range interactions 

Previous in vivo studies using DENV2 Renilla luciferase 
(Rluc) reporter replicon in BHK21 cells have shown that 
the 'flipping' mutation of PK1 severely reduced viral 
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Figure 4. Chemical probing of PK2 F i ip mutant RNA. The description of the scheme follows the convention of the Figure 1 legend. The flipping 
mutation within PK2 is marked by a blue rectangle. The insert in the left top corner represents increased Pb 2+ -induced hydrolysis within PK2 and 
TL1 regions. 



replication, whereas mutating both TL2 and PK1 to 
maintain their base pairing generated even more severe 
phenotype (8). Also, introducing compensatory mutations 
in the 5' CS (mutant TL2 F i ip PKl Flip -3) to restore base 
pairing between the 5'— 3' CS failed to improve replication, 
from which strong structure- and sequence-specific 
requirements for TL2 and PK1 function in DENV RNA 
replication have been proposed (8). Consequently, we 
sought to investigate the structure-function relation of 
both regions in the context of DENV-MINI RNA. 

The NMIA sensitivity profile of PK1 Piip mutant 
indicated extensive alternations within the 5'- and 3'- 
UTRs (Figure 5 and Supplementary Figures S7 and S8). 
The 5' UAR no longer formed a double-stranded region, 
but instead, it became a part of the novel SLB hairpin 
(G70-C105) with a reactive apical loop (G86-G89) and 
internal loops (U73-A77 and G99-U102). The 5' DAR 
region, previously double-stranded in DENV-MINI, 
rearranged into a reactive single-stranded loop (CI 05- 
G110). In addition, the PK1 'flipping' mutation abolished 
base pairing between 5' and 3' CS, disrupting the previous 
double-stranded 5'-3' CS. The 5' CS (C135-G144) now 
formed a reactive single-stranded region downstream of 
cHP, whereas the 3' CS (C613-U622) and the 'flipped' 
PK1 region (C613-A617) comprised a reactive loop down- 
stream of the 3' DB. In addition, the two hairpins previ- 
ously formed within the sequence of the 5'-end of the 
capsid protein gene were fused into an extended hairpin 
(G153-C197) with a weakly reactive internal loop (A158- 
G163 and A187-A192) and an A-C mismatch. 

Further alternations were noted within the 3'-UTR, 
where the 3' SL became a part of extended hairpin 
(A641-U719) proceeded by short hairpin (sHP) 



(C627-G640) containing moderately reactive GN«/RA 
hexaloop (Figure 5 and Supplementary Figures S7 and 
S8) (23,38). The sHP was shown previously to fold 
locally in the linear conformation of the DENV RNA, 
governing the balance between alternative DENV RNA 
structures (39). The 3' UAR, which no longer interacts 
with the 5' UAR, now formed both the rightmost side of 
the sHP and the leftmost side of the 3' SL (C638-U650), 
whereas the 3' DAR became a leftmost part of the sHP 
(C625-G630). As the 'flipping' PK1 mutation affects 
the secondary structure of neither the 3' nor the 5' DB, 
this result likely ruled out the existence of the previously 
proposed TL2/PK1 pseudoknot interaction (13,16). The 
VR region was not affected (C232-U448), with the excep- 
tion of a short helical region (A190-U194 and 
A443-U448) that in DENV-MINI RNA linked both 
DBs and VR domains, and now was predicted to be 
single-stranded and NMIA reactive. 

Pb 2+ ion-induced cleavage pattern for mutant PKlpup 
substantiated the SHAPE-predicted alterations of RNA 
secondary structure (Figure 5 insert and Supplementary 
Figure S9). Residues of SLB apical loop (A85-A88), as 
well as internal loop (U73-A78 and U99-U100), were 
readily hydrolyzed, whereas the A78-G84:U91- 
U98 stem was not cleaved. Also, the DAR region 
(A107-C109) and 5' CS (A137-C142) were susceptible to 
cleavage, which agreed with their single-stranded charac- 
ter. With regard to the 3'-UTR, moderate hydrolysis 
was noted within the stem of the extended 3' SL, as it 
contained numerous bulges and loops that relaxed its 
overall stability. Pb 2+ cleavage was most pronounced 
within the apex of the 3' SL (A674-C675). The 3' CS 
region containing the flipped PK1 exhibited 
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Figure 5. Chemical probing of PKl F ( ip mutant RNA. The description of the scheme follows the convention of Figure 1 legend. The flipping 
mutation within PK1 region is marked by blue rectangle. The insert in the left top corner represents increased Pb 2+ -induced hydrolysis within 
'flipped' PK1 residues and the novel motif, stem-loop B (SLB). 



single-cleavage points (A604, A605 and A620), 
whereas the apex of sHP (GN«/RA hexaloop) showed 
no susceptibility to hydrolysis. Overall, the PK1 'flipping' 
mutation disrupted the 5-3' terminal interactions by 



affecting the conformation of RNA motifs on both 
termini. 

Regarding the TL2 Flip PKl Flip mutant, the addition of 
'flipping' mutations within TL2 was designed to preserve 
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complementarity between TL2 and PK1. Also, SHAPE- 
assisted prediction of RNA folding showed that PK1 
was embedded in single-stranded region available for 
pairing with TL2. It was suggested previously that 
releasing PK1 from the 5—3' CS region provides the 
switch between the translation and replication mode of 
the DENV genome by rendering PK1 available for a 
pseudoknot interaction with TL2 (13). In our study, the 
RNA secondary structure predicted for TL2 F | ip PKl F |i p 
mutant resembled the structure of the PKl Flip mutant 
(Supplementary Figures S10 and Sll). Residues of 
neither PK1 nor TL2 displayed a change in NMIA sensi- 
tivity, contradicting the implication that TL2 is involved 
in long-range interactions with PK1 in the context of 
DENV-MINI construct. Pb 2+ -induced hydrolysis of 
TL2 Flip PKl F | ip mutant provided a pattern of cleavage 
comparable with that obtained for the PKl F i ip mutant 
(Supplementary Figure SI 2), the only difference being 
lower intensity of hydrolysis within the internal loops 
and rightmost apical loops of both DBs. 

We also investigated how introducing 3-nt compensa- 
tory mutations within the 5' CS (U141-C144) to restore 
the 5—3' CS in the presence of 'flipped' TL2 and PK1 



affected DENV MINI RNA structure. SHAPE analysis 
of mutant TL2 F i ip PKl F i ip -3 predicted almost complete 
reversion to the wild-type secondary structure (Figure 6 
and Supplementary Figures SI 3 and SI 4). The only re- 
arrangement occurred within the 5'-3' CS region, which 
was shortened from 11 -nt helical region to 6-nt stem 
(C135-U140 and A624-A618), releasing PK1 (C613- 
A617) and compensatory residues into single-stranded 
loops. Thus, the complementarity and single-strand char- 
acter of both the downstream component of 5' CS and 
PK1 proved insufficient to maintain their base pairing 
interaction. Pb 2+ -induced hydrolysis of TL2 Flip PKl Flip -3 
mutant (Figure 6 insert and Supplementary Figure 15) 
indicated moderate cleavage within the region carrying 
compensatory mutations, as well as within PK1 
(A615-C616). 

Collectively, chemical probing analysis of DENV-MINI 
mutants supported the notion that, rather then being 
involved in the TL2/PK1 pseudoknot interaction, the 
structure- and sequence-specific composition of the PK1 
region is essential for sustaining the correct conformation 
of the 5' and 3' terminal motifs necessary for 5—3' terminal 
RNA-RNA interaction. 




Figure 6. Chemical probing of TL2 F ii p PKl F ii p -3mutant RNA. The description of the scheme follows the convention of Figure 1 legend. The 
3-nt compensatory mutations introduced within 5' CS, the 'flipping' mutations within PK1 and TL2 region are marked by blue rectangles 
and trapezoid. The insert in the left top corner represents the pattern of increased Pb 2+ -induced hydrolysis within 'flipped' PK1 and compensatory 
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3D modeling of DENV-MINI TL1/PK2 pseudoknot 
interactions 

A 3D structural model of the DENV-MINI RNA was 
generated using RNAComposer, a fully automated 
server that predicts 3D structural models from RNA sec- 
ondary structure (40) (Figure 7). As the server does not 
accept sequences longer than 500 nt, a 474-nt derivative 
sequence was created by deleting the VR domain and 
closing the remaining short helical region (A190-U194 
and A443-U448) with a -GAGA- tetraloop. A dot- 
bracket notation generated by the RNAstructure 
software was adjusted to account for this internal 
deletion, as well as the predicted TL1/PK2 pseudoknot, 
and subsequently provided to RNAComposer. Ten 3D 
RNA models were generated and analyzed, taking into 



account their secondary structure topology, sequence 
homology, structure resolution and free energy. In 
addition, the quality of predicted models was evaluated 
using the MolProbity tool (24). The best of these models 
was manually adjusted to correct for an unlikely fold 
within PK1, and the affected regions were subjected to 
energy minimization (Accelrys Discovery Studio v3.5) to 
create the structure of Figure 7. 

The model shows close association of the 5'- and 
3'-UTRs (Figure 7A and B), in agreement with previous 
studies, suggesting that efficient transfer of RNA-depend- 
ent RNA-polymerase from the 5' SLA to the 3' SL motif 
is essential for initiation of (— ) strand RNA synthesis (12, 
17). In addition, positioning of the hybridized 5—3' 
termini and the TL1/PK2 pseudoknot on opposite 
termini of the DENV-MINI RNA is consistent with the 




Figure 7. Proposed 3D model of truncated DENV-MINI RNA. The internally deleted DENV-MINI RNA is depicted (A) in its entirety or 
(B) highlighting the proximal 5'- and 3'-UTRs. Specific RNA domains are color-coded as shown in the key. The cis-acting motifs are labeled as 
follows: SLA TL, stem-loop A top loop; SLA SSL, stem-loop A side stem-loop; 3' SL, 3' stem-loop; dsUAR, double-stranded upstream of AUG 
region; PK, pseudoknot; TL, top loop. For clarity, the remainder of the RNA molecule is colored in gray. 
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two domains playing independent structural/regulatory 
roles in DENV replication. TL2 of the 3' DB is likewise 
spatially separated from other structural motifs, suggest- 
ing a functional role for TL2 that does not involve inter- 
actions with other RNA elements in the context. The 3' 
DB itself is linked to adjacent segments of the DENV- 
MINI RNA by relatively long, presumably flexible 
segments of single-stranded RNA. It is conceivable that 
a second TL2/PK1 pseudoknot could be formed by 
rotation of the 3' DB motif, thereby bringing TL2 into 
proximity with PK1. However, such a rearrangement 
would also require partial dissolution of the hybridized 
5—3' termini. It is important to note, however, that there 
is no evidence that base pairing between TL2 and PK1 
occurs in the DENV-MINI RNA in vitro. In general, 
although caution must be exercised in interpreting these 
3D structural data, the model presented here serves to 
provide plausible insights into how RNA substructures 
may interact during DENV replication. 

DISCUSSION 

Following infection, the (+)-stranded DENV RNA must 
serve as a template for replication, an mRNA for transla- 
tion and a substrate for encapsidation. The function(s) 
assumed by these RNAs are dictated in part by the 
specific local and long-range RNA interactions investigated 
in the present study. Our combined SHAPE, aiSHAPE and 
Pb 2+ -induced cleavage analysis of the DENV-MINI RNA 
establish the existence of multiple motifs required for initi- 
ation of (— ) strand RNA synthesis (Figure IB and 
Supplementary Figure SI). The presence of the 5' SLA, 
demonstrated crucial for initial binding of RdRp (6,17), 
the single-stranded polyU tract providing the flexibility to 
the 5'-UTR (17), together with the rearrangements 
occurring within 3' SL (17,30), established replication 
competency of DENV-MINI RNA (Figure IB and 
Supplementary Figure SI). Our results have experimentally 
verified long-range interactions between the 5' and 3' 
terminal regions, which involved hybridization between 
(i) 5—3' CS, previously demonstrated to initiate the 5—3' 
long range RNA-RNA hybridization (5,15), (ii) the 5'-3' 
UAR, required for the correct orientation of the SLA with 
respect to the 3' end (30) and (hi) 5'-3' DAR, extending the 
CS-initiated 5—3' interaction (11). We have predicted add- 
itional regulatory motifs that might contribute to the 
overall functionality of DENV-MINI RNA. Probing 
analysis identified tandem DBs spaced by A-rich single- 
stranded regions in the 3'-UTR, as well as internal loops 
and hairpins with embedded GNRA/GNRA-like motifs 
within the VR domain (Figure IB and Supplementary 
Figure SI). These elements may represent candidate 
assembly sites for proteins, as previous studies identified 
several host factors that not only bind to DENV UTRs 
but also actively participate in regulation of replication 
and translation efficiency (41^14). In addition, -GNRA- 
tetraloops could act as 'anchors' by interacting in a 
sequence-specific manner with motifs at distant locations 
along the molecule to enhance RNA self-folding (45,46). 

We have provided here novel insight into tertiary inter- 
actions that exist within the DB structures. Previous 



phylogenetic comparison of four sero-groups of flavi- 
viruses suggested formation of two pseudoknots, TL2/ 
PK1 and TL1/PK2, involving base pairing between the 
DBs TLs and the downstream complementary PK 
regions (13). These authors suggested that the single- 
base pair difference (A-G), distinguishing PK2 and 
PK1, conserved in all DENV strains, must have a func- 
tional reason to differentiate between both motifs (13). 
Manzano et al. (8) also indicated dissimilar functional 
roles of TL2/PK1 and TL1/PK2 in DENV replication 
and translation. 

Our SHAPE-directed probing implied that the base 
pairing interaction between TL1 and PK2 formed the 
core of an extended, H-type pseudoknot, as both 
regions, even though embedded within single-stranded 
loops, contained unreactive residues (Figures IB and 2). 
Also, both sequences were resistant to the complementary 
strategy of Pb 2+ -induced hydrolysis, supporting their 
potential base pairing interactions (Figure IB insert and 
Supplementary Figure S2). aiSHAPE lends further 
credence to TL1/PK2 tertiary interactions (Figure 3 and 
Supplementary Figure S3), as outcompeting one base 
pairing partner increased reactivity of the other and vice 
versa. Moreover, SHAPE analysis of the PK2 F | ip mutant 
indicated a major increase in NMIA reactivity of TL1 and 
PK2 nucleotides, validating disruption of the TL1/PK2 
base pairing interaction (Figure 4 and Supplementary 
Figures S4 and S5). Consequently, we propose that in 
the context of DENV-MINI RNA, the TL1 loop of the 
5' DB is involved in hairpin-type (H-type) pseudoknot 
interactions with the downstream PK2 region (Figure 2) 
(47). This pseudoknot was predicted to have two helical 
regions, including base paring between PK2 (G526-C530) 
and TL1 residues (G470-U474), as well as the stem of the 
leftmost hairpin (C462-G466 and C478-G482), and three 
single-stranded loops: A468-G470, U474-C476 and 
C516-C525 (Figure 2B). 

We have also provided a 3D model of the TL1/PK2 
interaction using the recently reported program 
RNAComposer (40). As noted in Figure 7, the hybridized 
5'-3' termini and the TL1/PK2 pseudoknot are located on 
opposite ends of the DENV-MINI RNA, suggesting an 
independent structural/regulatory character. Also, the 
fact that hybridization of LNA/DNA chimeras to either 
TL1 or PK2 failed to induce extensive structural changes 
within the interacting 5'-3' termini (Supplementary Figure 
S3) suggests the TL1/PK2 pseudoknot might indeed 
function as an autonomous element. In addition, TL2 of 
the 3' DB is spatially separated from other structural 
domains and linked to adjacent segments of the DENV- 
MINI RNA by relatively long single-stranded regions. It 
is possible that a second TL2/PK1 pseudoknot could be 
formed by rotation of the 3' DB, thereby bringing TL2 
into proximity with PK1. However, such a rearrangement 
would also entail partial dissolution of the hybridized 5'-3' 
termini and perhaps even the TL1/PK2 interaction, such 
as might occur in the presence of viral proteins and/or 
host factors, or in the event the 5' terminal portion of 
the genome has been removed by ribonucleases (10) or 
is sequestered during translation. A similar hypothesis 
was proposed based on MPGAfold-directed secondary 



Nucleic Acids Research, 2013, Vol. 41, No. 9 5087 



structure prediction of DENV MINI RNA, which 
indicated differences in the folding behavior of the two 
DBs (8). In the lack of the 5-3' communication, it was 
proposed that TL2 could be available to bind PK1 in the 
3' DB to form the putative TL2/PK1 pseudoknot (8). 

Pseudoknots are associated with a remarkably diverse 
range of biological activities, including translational 
control via internal ribosome entry site (48,49), ribosomal 
frameshifting (50,51), tRNA mimicry (47,52-54) and 
temporal control of the termination of translation and 
initiation of genome replication (55). It has been shown 
that 'flipping' of PK2 in the context of luciferase-based 
DENV2 replicon did not affect translation, but disrupted 
replication (8), from which we conclude that the TL1/PK2 
pseudoknot participates in regulation of DENV RNA 
synthesis. It is also possible that the TL1/PK2 pseudoknot 
might create a critical recognition element for host 
proteins. Supporting this notion, pseudoknots have been 
shown to cause distortion of A-form helices, creating a 
bend in the complex to provide a recognition site for 
ligands (56). Also, it was shown that DEAD-box RNA 
helicase (DDX6) binds to both DBs and their correspond- 
ing PKs regions (43). These authors suggested that the 
putative pseudoknots might provide assembly sites for 
DDX6, and that proper assembly of this protein had func- 
tional consequences for DENV replication. 

Our combined probing and mutational analysis also 
addressed the role of PK1 in the DENV-MINI RNA sec- 
ondary structure. The 'flipping 1 mutation introduced into 
PK1 significantly rearranged the conformation of the 
5'- and 3'-UTRs (Figure 5 and Supplementary Figures 
S7, S8 and S9). The 5' and 3' CS became single- 
stranded, whereas the UAR and DAR, previously 
required for 5—3' communication, became involved in for- 
mation of the 5' SLB, the 3' sHP and an extended 3' SL 
motif. The presence of 5' SLB was proposed as specific for 
the linear conformation of DENV RNA (39), reducing the 
ability of the 5' UAR to hybridize to the 3' partner (12). In 
addition, the sHP was suggested to fold locally in the 
linear conformation of DENV RNA (39). Although its 
function is uncertain, the sequences within the sHP were 
proposed to regulate the equilibrium between the DENV 
RNA alternative conformations for organizing multiple 
functions of the viral genome. 

SHAPE-directed prediction of RNA folding indicated 
similar structure rearrangements for the TL2 F ij p PKl F ii p 
mutant RNA (Supplementary Figures S10-S12), suggest- 
ing PK1 is mostly responsible for correct long-range com- 
munication between the 5' and 3' terminal regions. Even 
the presence of compensatory mutations within 5' CS 
(TL2 Fli pPKl F i ip -3 mutant) to restore base pairing with 
PK1 (Figure 6 and Supplementary Figures S13-S15) did 
not entirely reinstate the secondary structure of DENV- 
MINI RNA. Previous studies using DENV-MINI RNA 
as templates in in vitro RdRp assays have shown that rep- 
lication of PKl Flip , TL2 Flip PKl Flip mutants was signifi- 
cantly affected, supporting the importance of the 5—3' 
CS interaction for RNA synthesis (8). However, the intro- 
duction of compensatory mutations at the 5' CS in 
DENV-MINI TL2 Flip PKl Flip -3 mutant restored replica- 
tion to nearly wild-type level. This was in contrast to the 



effects of the PK1 'flipping' mutations on replication of 
luciferase-based DENV2 in BHK-21 cells, which showed 
essentially no replication (8). Related studies using 
luciferase containing Dengue virus (15), Kunjin virus 
(57) and West Nile virus replicons (58) carrying point mu- 
tations at the 3' CS are consistent with our observations. 
In all cases, the altered sequence of 3' CS attenuated rep- 
lication, even when mutations were compensated at the 5' 
CS. Also, a mutation introduced into the 3' CS extending 
the region of complementarity between 5-3' CS by one 
nucleotide was shown to negatively affect replication 
kinetics of the subgenomic DENV reporter replicon (59). 
Apparently, not only the secondary structure of 5—3' CS 
region but also its sequence composition governs the 
efficiency of RNA synthesis. 

In addition, the TL2 F | ip PKl F i ip mutant (Supplementary 
Figure S10) provided insight into the second, previously 
proposed TL2/PK1 pseudoknot (13,16). It was suggested 
previously that PK1 could modulate DENV genome cyc- 
lization via either a 5—3' CS interaction or by base pairing 
with TL2 (13,25). Also, studies showed that the long- 
distance 5'-3' CS complementarity interferes with efficient 
translation of DENV2 luciferase (LUC) reporter system in 
Vero cells (25). It was suggested that the 5'— 3' CS inter- 
action may play role in regulating the transition from 
the translation to the replication phase of the infection. 
Thus, we expected that release of PK1 from 5—3' CS 
together with the conserved complementarity between 
both partners, would support formation of TL2/PK1 
pseudoknot. SHAPE analysis, however, indicated no 
reactivity change within 'flipped' TL2 and PK1 
(Supplementary Figures S10 and Sll). As the addition of 
'flipping' mutations within TL2 did not cause any further 
structural rearrangements compared with PKl Flip mutant 
(compare Figure 5 and Supplementary Figures S7 and S8), 
this would dismiss the involvement of TL2 in tertiary inter- 
actions. However, we cannot exclude the possibility that 
in vivo a TL2/PK1 interaction could be stabilized by par- 
ticipation of host or viral factors for the purpose of other 
molecular processes, i.e. translation, before the 5'-3' CS 
interaction. Indeed, it has been proposed that the viral 
NS3 helicase can modulate the folding of ew-acting RNA 
elements via dual RNA unwinding/annealing activities, 
providing the switch between alternative conformations 
of the DENV genome (38,60). 

Considering that DENV RNA must assume alternative 
conformations (39) to fulfill a continuous need to function 
as mRNA, the replication template, as well as the viral 
genome, we must bear in mind that the structure 
reported here may be only partially representative of the 
in vivo viral RNA population. Further studies will address 
folding of DENV RNA at different stage of viral life 
cycle to understand the mechanisms that modulate multi- 
functionality of viral RNA and, as a result, promote 
rational design of live attenuated vaccines and/or antiviral 
therapeutics. 
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